Isotropic dispersion, line shape, and remnant Fermi surface in one hole problem 
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It is shown that the isotropic dispersion relation, line shape of the spectral function, and remnant 
Fermi surface found in recent photoemission experiment in insulating Ca2Cu02Cl2 by Ronning et 
al. (Science, 282, 2067 (1998)) can be consistently explained by a direct holon hopping process 
which fundamentally undermines the conventional self-consistent Born approximation approach to 
the t — J model. 

74.20.Mn, 79.60-i, 71.27.+a 



The problem of one hole moving in the antiferromag- 

[ netic (AF) spin background can give us some impor- 
tant clue about the finite doping. For many years, the 

] self-consistent Born approximation (SCBA) method 0-^] 
based on the t — J model has provided a standard de- 

] scription of the one hole problem, consistent with exact 
numerical calculations j|] up to 32 sites which pre- 

\ diets a sharp quasiparticle coherent peak in the spectral 
function with an anisotropic dispersion relation around 

] the Fermi points ko — (±7r/2, ±7r/2). 

But the SCBA results are not consistent with angle- 
resolved photoemission spectroscopy (ARPES) measure- 
ments in the parent compounds of high- Tc superconduc- 

; tors: SrsCuOaCla @,§ and CazCuOaCla §, which in- 
dicate that while the "quasiparticle" peak is not sharp 

' at all with an intrinsic broad energy structure, its dis- 
persion is fairly isotropic around the Fermi point kg. A 
further challenge to theory comes from the recent high- 
quality measurement |^ revealing the so-called remnant 

^ Fermi surface (RFS) in the electron relative momentum 
distribution rtj^ , in contrast to the low-energy Fermi point 
structure. 

To explain the isotropy of the dispersion relation, the 
second and third neighbor hopping t' and t" adding to 
the t — J model were suggested in literature |p|, p^ , pT[ 
which could improve the agreement between the theory 
and experiment under some proper choice of parameters. 
But Laughlin has proposed a different view in which 
the isotropic dispersion may come from a more intrinsic 
origin: the spinon spectrum. In this spin-charge separa- 
tion picture, the size of the dispersion is also accounted 
for by the superexchange J without t' and t" , and the in- 
jected hole decaying into a spinon-holon pair can provide 
' a natural broad structure in the spectral function. Such 
a spin-charge separation scenario is fundamentally differ- 
ent from the Landau quasiparticle picture in the SCBA 
approach. 

The purpose of this paper is to reexamine the one hole 
issue based on the t — J model. We show that due 
to the existence of a direct holon hopping term miss- 
ing in the SCBA scheme, the low-energy physics can 
be completely changed from the SCBA prediction, and 



the isotropic "quasiparticle" dispersion, line shape of the 
spectral function, as well as the remnant Fermi surface 
structure observed in experiment can all be explained as 
the consequence of the spin-charge separation, which in 
some main aspect is consistent with the proposal made 
in Ref. |l|]. 

We start with the slave-fermion representation da = 
f}ha §. In the t - J model Ht-j = Ht + Hj, 
Hj = -i Efe)(Af,)^Af^. with Af^. = E [Here 
bia- = {—aybia, with (— <t)* explicitly keeping track of 
the Marshall sign.] The motion of a doped hole is gov- 
erned by Ht = -tY,(ij)Bjif}fj + H.c. with 13 = 

cr^ja^icr- Evcu though the Schwinger-boson mean- 
field theory (SBMFT) ^ works quite well at half- 
filling, characterized by the bosonic resonating-valence- 
bond (RVB) order parameter {^fj), the one hole problem 
is highly nontrivial. The reason can be attributed to the 
fact that the hopping integral (Bji) = in Ht- In the 
SCBA approach, the hole has to acquire its kinetic en- 
ergy through the dynamic fluctuations of Bji as follows. 
To the leading order in a large-S (spin wave) expansion, 
one may write Q 

B.^^B^f"^ ^boJ2cT[bl + -b,a], (1) 

where Bq denotes the condensed part of the Schwinger 
boson field, corresponding to AF long-range order 
(AFLRO) Then a hole created by Ci^ = //&o(-cr)* 
can propagate by an SCBA assisted hopping and behave 
just like a Landau quasiparticle with four Fermi points 
atko §f . 

However, an important process is missing in the above 
picture in the spin-1/2 case. We will show below that 
due to the RVB nature of the SBMFT description 
a direct hopping term for the holon actually exists but 
is omitted in (p. For the sake of clarity, one may set 
the temperature T at O"*" where bo (and the AFLRO) 
is switched off such that (|l|) vanishes. Then one can 
rewrite Bji in the following form Bji = {Bj^) e**J' where 

B% = EJl'h. and 4>r^ = ±i^/m " ^n] = ±1 
depending on ]" or | spinon exchanged with the holon at 
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the given link (ij)]. If (B^i) = 7^ 0, the following bare 
holon hopping term will emerge 



(2) 



with teff = tB^. Here (Bji) = is entirely due to 
^e^J^ji) = but the correct measuring of the gauge phase 
strength should be the gauge-invariant quantity e*Sr 
for a closed path F. One may write e*^r'^^' = (— l)^r 
with denoting the total number of | spins (or t spins 
by symmetry) "exchanged" with the hole along the path 
r. Here (— l)^r is nothing but the phase string factor 
first discussed in Ref. in a rigorous formulation. 

The proof of i?° 7^ is straightforward. At half-filling, 
one has the Bogoliubov transformation 



Uk7k<T - -fkTk-CT 



(3) 



in which 7kcr represents an SBMFT elementary excitation 
operator, with a gapless spectrum at T = 0+ Q 



Ek — 2.32J^ 



(4) 



and u^. — U] 



1, 2MkWk = Ck/Vl ~ ^k where = 
(cos kx + cos ky)/2. It is crucial to recognize that the fac- 
tor TykCT in (H) which satisfies r]a = rf_^ and \rj\ = 1 cannot 
be completely determined at the mean-field level. Specif- 
ically, the mean-field order parameter {bi„bj-a) is inde- 
pendent of different choices of rjka ■ Then one may exploit 
this symmetry for the one hole case. At each fixed hole 
position, there is a sub-Hilbert-space for spins which at 
the mean-field level is described by the Bogoliubov trans- 
formation (^. By choosing Tyko- = [~S(7n(^k)]'^'' where 
kh = if the hole is on the even sublattice and kh = 1 
if the hole is on the odd sublattice site, it is easy to 
find that = Sk^o l^kl^k ~ 0.4 in Ht where a hole 
changing sublattices is involved, while the mean-filed re- 
sult remains the same for the superexchange term Hj. 
Note that such a hole-position-dependent choice of the 
phase 77kCT is perfectly legitimate as in the sub-Hilbert- 
space in which the hole position is given, the commuta- 
tion relations of b, b^ remain unchanged. 

The existence of a bare hopping term Hh in (j2|) will 
fundamentally change the SCBA result. Note that one 
may express Qo- = + ha with ha = // : ha ■ 
with : ■= bier — &o(— c)*- Due to (|^), the propaga- 
tor for Cio- has a leading contribution Gg ~ iGf ■ Gb 
which refiects the spin-charge separation, where Gf ~ 

-m{fKt)fM)) and Gb = -i{Tt{: ha{t) :: 5],(0) :)). 
On the other hand, the Hamiltonian corresponding to 
the linear spin- wave approximation (fil) may be rewritten 



as HJ^^^ — tJ2{ij)a (clffCjtT 



+ 5,1 c 



H.c, which 



literally shows a process for the quasiparticle Ca to de- 
cay into a holon-spinon pair represented by Co-. When 
both holon and quasiparticle have no bare kinetic energy, 
HJ^^^ represents a virtual process and an SCBA treat- 
ment of it reproduces the well-known results • But 
once the holon acquires a bare kinetic energy whereas 
the quasiparticle has none, H^^^ simply describes the 
break-up process of a quasiparticle, where the propaga- 
tor Ge for Cct should become proportional to Gg to the 
leading order approximation. Note that any higher order 
corrections to the quasiparticle self-energy also add to the 
holon self-energy so that the initial instability is always 
maintained. Therefore, the SCBA approach is totally 
undermined once the spin-charge separation is present 
where the holon as a free particle has its bare kinetic 
energy. 

The spin-charge separation reflected in Ge leads to the 
convolution law of the spectral function in the following 

form 



(5) 



where the r.h.s. is obtained after using the SBMFT re- 
sult, ImGh = —'k[u^S{uj — Ek) — v^^d^uj + -Ek)], for the 
spinon part, while the holon spectral function ^'*(k, oj) ^ 
at > which is to be determined by Hh in (||). By 
noting that 



giko-(rj-ri)g- 



(6) 



with S4>ji = ±(7r/2)crjj and {S(f>ji) = 0, one finds that 
yl''(k, w) should be peaked at k = ko at low energy, and 
correspondingly the spinon spectrum £^k+ko will essen- 
tially determine the low-energy "edge" of the electron 



spectral function in 



Here we note that 



will 



generally play a role causing the incoherent motion of 
the holon as opposed to the coherent spinon with a i5- 
function-like spectral function. As an approximation, we 



first simply treat 



as a random fiux field and de- 



termine A'^CkjOj) numerically. We note that the overall 
shape of ^"^(k, uj) will not be sensitive to such an approx- 
imation as to be discussed below. 

The typical behavior of A'^{'k,Lu) versus uj is shown in 
Fig. 1 at different k's along (0, 0) to ko. The small arrows 
mark the energy of i?k+ko which is indeed correlated with 
the position of the "coherent" peak or more precisely the 
edge of A^(k, oj). The line shape of the spectral function 
in Fig. 1 reflects essentially the convolution law of spinon 
and holon in a spin-charge separation context, which is 
strikingly similar to those found in the ARPES measure- 
ments In particular, the spinon spectrum -Ek+ko 
fits the ARPES data reasonably well over the whole Bril- 
loiun zone as shown in Fig. 2 along kg (S) to (0,0) (F) 
and (0,7r) (X) (which are symmetric in our theory) with 
J = O.lAeV. A very flat (only about 13% change) fi'k+ko 
along (0,0) to (0, tt) is also in good agreement with the 
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experiment. Therefore, the theory can naturally account 
for the ARPES without t' , t" terms. 

In Fig. 1, the holon spectral function is computed at 
teff = 2J « 5 J) with the random (white noise) flux 
strength chosen between ±0.l7r per plaquette in 48 x 48 
lattice. But we emphasize that the general line shape in 
Fig. 1 looks insensitive to those parameters which mainly 
affect A'^. This is because the broad feature of the spec- 
tral function is chiefly due to the convolution law which 
reflects the spin-charge separation and the high en- 
ergy part of A"^ is essentially decided by the bandwidth 
Steff) of the holon spectrum, not the sharpness of A^ 
itself. Note that the coherent factor v'^ and the "incoher- 
ence" of the holon caused by the gauge phase in (||) are 
crucial for the coherent spinon spectrum to "show up" in 
A'' (k, oj) . In the one-dimensional case, the phase string 
also plays a role to shape the one-hole spectral function 
as discussed in Ref. [|l6| . 

Beyond the random-flux treatment, the precise form of 
S4>ji will be important in determining the relative change 
at different momentum of an energy- integrated A'^: the 
so-called relative momentum distribution |^ deflned by 
n^([jo) — ^^'^{k,uj). Here if uiq is taken to — oo, 
will reduce to the momentum distribution n^. which 
should be featureless (= 1/2) for the one- hole case [note 
that (cl^Cja) = (bl^fifjbja) = Sij{bl^bia) using the non- 
double-occupancy constraint]. Thus, within the frame- 
work of the t — J model, the observed RFS feature in 
^ for a flnite luq (~ 0.5eV) has nothing to do with the 
half-flUing ground-state and must solely come from the 
dynamics of the doped hole. We have already seen that 
the characteristic "Fermi-point" ko arises from the lead- 
ing term of the gauge phase in and we will see 
that beyond the random flux approximation S(j)ji may be 
responsible for an additional momentum structure in n^. 

The gauge phase (f> can be more precisely described as 
a nonlocal phase string effect as noted before. It may be 
clearly seen by formulating G/ according to Ref. |]l^ as 
follows 



Gf{J,^■,E)= (5,(K};£;))(-1)^^ 

{c}{<t4 



(7) 



with Sr — 



Ut,Gj{t) 



6t 



Gj, Gj — 



1/{E - Hj), and = ^^(1 - as)/2. In S'c, the site 
nis is at a path c that the holon travels, which connects 
sites i and j: rrio = h ttli, "m-K = j- In particular, one 
can prove ||l^ that {Sc{E)) > at low energy such that 
the phase string effect introduced by (—1)''^= is basically 
nonrepairable. What we need to know is how such an os- 
cillating phase string shapes the high energy part of A"^ to 
determine n^. Due to its nonrepairable nature, the sin- 
gular part of the phase string effect can be explicitly built 
into the wave function through a unitary transformation 



the elementary charge and spin excitations whose details 
are not important to our problem below. The key thing 
is the phase shift Q'^^""^ which will then modulate the 
electron propagator by 



1 = p"* 



with A-'' (r) 



zx(r-r,) 



la, leading to a, 



X e 



Here ain- refers to 



(8) 

in the equal- 
time limit. Such a phase string factor originated from (pji 
can be similarly evaluated using the method outlined in 
Ref. ijl^ , which leads to four "incommensurate peaks" at 
(±7rK, 0) and (0, zLttk) in momentum space with k ~ 1 
||l8| . If this oscillating factor solely decides the momen- 
tum structure of nj^ with the rest term in the propagator 
mainly contributing a broadening, then one gets a con- 
tour plot of in Fig. 3 in terms of the superposition 
of aforementioned four peaks (with k = 0.75) at an ar- 
bitrary broadening (which is presumably controlled by 
the energy cutoff luq: when loq — > — cxd, the broadening 
in the momentum space should also go to infinity such 
that rik — 1/2). Note that the values of here are only 
for relative comparison. We see that the overall feature 
is in qualitative agreement with the experimental data 
Q whose one quarter portion is also shown in the inset 
of Fig. 3 for comparison. Experimentally |^ the RFS is 
defined at the sharp drop of n^, in k-space. 

Finally, we point out that strictly speaking the AFLRO 
and RVB condensations only exist in the thermodynamic 
limit and both (^ and (||) do not directly apply to a 
finite system. In a small sample, both amplitude and 
phase fluctuations in Bji are expected to be important 
and cannot be distinguished like in (|^) where 5° 0. 
And in this case a generalized SCBA scheme may be con- 
structed to account for the leading order effect in the hole 
dynamics. This should be the reason why the small-size 
exact calculations have always identified a quasiparticle, 
consistent with the SCBA as well as power-method re- 
sults 1^,^,^. However, such an SCBA quasiparticle is 
expected to eventually break up into holon-spinon com- 
posite at large lattice so long as there is a bosonic RVB 
condensation in the thermodynamic limit to warrant 

(I). 

In conclusion, a direct holon hopping term is identified 
in the one hole problem due to the RVB spin pairing in 
the SBMFT framework, which is omitted in the SCBA 
approach. Consequently a drastically different character- 
ization of low-energy physics emerges: The SCBA quasi- 
particle becomes unstable and is replaced by a holon- 
spinon composite in which the holon is incoherent due 
to the phase string effect associated with the bare hop- 
ping, while the coherent spinon spectrum determines the 
isotropic "quasiparticle" dispersion relation. The line 
shape of the spectral function simply reflects the spin- 
charge decomposition, while the remnant Fermi surface 
is caused by the oscillating phase string effect at high 
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energy. These properties well explain the recent ARPES 
measurements within the framework of the t — J model. 
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Fig. 2. The quasiparticle spectra determined by the 
ARPES 1^ (open square) and by the present theory (solid 
curve). 

Fig. 3. The contour plot of the electron momentum 
distribution exhibiting a "remnant Fermi surface" struc- 
ture. The experimental data ||^ are shown in the insert 
at the upper right corner. 
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